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Introduction 


Ce projet vise à étudier une installation électrique domestique existante. Le choix 
de l'appartement sur deux étages dans lequel je vis actuellement s'est imposé 
comme plus pratique et comportant suffisamment d'éléments pour avoir le 
mérite d'être étudié. 


Ce projet permettra d'aborder les notions de cours comme la sécurité électrique, 
le calcul des puissances actives, réactives, apparentes, des courants et tensions 
efficaces. De plus, des méthodes des réseaux de Kirchhoff et le théorème de 
Boucherot et de Tellegen seront utilisés pour les calculs, couplés à des méthodes 
de résolution numériques sous python. 


On commencera par décrire de manière détaillée l'installation électrique. On 
décrira ensuite certains éléments sécuritaires, reliés à la norme NF C 15-100 puis 
on réalisera les calculs des circuits électriques. Finalement, avant de conclure, on 
comparera la consommation électrique trouvée avec la consommation réelle. 


Afin de faciliter la clarté de la lecture, certains documents sont directement 
intégrés au texte et d'autres données en annexe. 


)) Description de l'installation 


A) De la production à la consommation d'électricité 


Le système étudié est localisé en rouge sur le schéma ci-dessous. Il s'agit d'une 
maison individuelle, installation en monophasée. L'électricité, avant d'arriver dans 


l'habitat, suit le schéma suivant : 


RÉSEAU DE TRANSPORT RÉSEAU DE DISTRIBUTION 
Grand transport Répartition Distribution 
Lignes haute tension Lignes moyenne tension Lignes haute 
et basse tension 
U > 50 kV 1 kV < U < 50 kV 50 kV < U < 1 kV 


Hôpitaux 
ps 
CELL Fuss 
éseau . 
SC X À. ferroviaire r\ CR 
Æ Eu" 
À + Æ 
/Æ\ "À À 
un 
CREER 
oral — "À — HE we 
Postes électriques 
(transformateurs) 
Système étudié 
Usines Zones 
de consommation 
Lieux de production Infrastructures (logements, 
spéciales commerces, etc.) 


(centrales électriques) 


Fig 1 : Production et transport de l'énergie électrique 


1) La production d'électricité, non représentée sur le schéma, est 


généralement assurée par EDF (même si depuis l'ouverture du marché de 


l'électricité à la concurrence en 2007, il existe de plus en plus d'autres 


acteurs décentralisés de production d'énergie) 


2) Le transport haute tension HTB (afin d'éviter les pertes par effets joules 
au maximum) est assuré premièrement par le réseau de transport 


électrique (Rte!) entre le lieu de production et les postes sources. Ces 


postes sources, par l'intermédiaire de transformateurs, convertissent 
alors la très haute tension HTB (U > 50 KV) en moyenne tension HTA (1 kV 


< U < 50 kV). 


3) Le transport haute tension HTA et basse tension BT est ensuite assuré 
par Enedis, filiale d'EDF se chargeant de distribuer la HTA aux 
infrastructures la nécessitant (Industries, Hôpitaux) et la BT aux 


particuliers (cas du système étudié). 


B) Eléments de facturation 


Grâce à la facture du fournisseur EDF, avec lequel un 
contrat est souscrit, on peut accéder à la puissance 
souscrite. On accède également à la puissance souscrite 
par le compteur Linky disposé à côté du compteur 
électrique. 


La puissance souscrite est de 6kVA. On en déduit 


immédiatement : 
S 6000 
Leff = Var = 30 — 26,09 A 
Cette valeur correspond donc au courant maximum 


soutiré par les équipements électriques domestiques. 


Fig 2 : Puissance souscrite de 6 


kVA affichée sur compteur Linky 


1 Le site Web de RTE a mis en place un outil de suivi “éCO2mix" très intéressant permettant de 
suivre en temps réel la consommation et la production d'électricité en France. Le site est 


disponible à l'adresse suivante : https://www.rte-france.com/eco2mix 


C) Description schématique de l'installation 


Afin de décrire l'installation, AutoCAD Electrical semblait particulièrement adapté 
en tant que logiciel le plus utilisé pour représenter des schémas électriques. 
Cependant, suite à l'exigence du moteur de la base de données Microsoft Access 
2016 pour faire fonctionner le logiciel, difficile à télécharger sur le Web, 
l'installation n'a pas pu être effectuée. Pour cette raison, la version standard 
d'AutoCAD a été utilisée. Pour davantage de clarté, deux plans ont été effectués, 
représentant les 12 pièces constitutives de la maison: 


1) Un premier plan architectural, représentant les différentes pièces de la 
maison et quelques équipements non-électriques. 


2) La légende est ensuite donnée, associée aux plans électriques simplifiés 


des deux étages, sur lesquels sont représentés uniquement les appareils 
fixes usuels les plus consommateurs d'électricité. 


Etage 1 Etage 2 


Chambre 2 Chambre 1 Chambre 4 


111 / 


Chambre 3 Chambre 4 
Couloir 1 


Couloir 2 


Cuisine 


Salon 


1 


vers étage 1 vers étage 2 


Le vues (E) 


vers étage 2 vers étage 1 ||*| 


Palier 


N\ Entrée 


Fig 3 : Plan architectural de la maison 


Le choix a été fait de décrire les appareils uniquement dans la partie 111) Analyse 
de la consommation électrique. L'ensemble des appareils est cependant 
également décrite dans le tableau général en annexe 2. 


Légende 
Element Description Nombre 
Prise de courant 

=} 28 

avec terre 

2 Double prise de 
:1] courant avec terre 5 
—TV Prise TV 2 
Tr Prise téléphone 3 
Interrupteur simple 10 
"il allumage 
F4 Interrupteur Va et 4 
Vient 
Fil électrique N+P+T 13 
Interrupteur-lampe 
() Lampe fluocompacte 3 
2 
a Lampe LED 11 
Appareil électrique à 

Ï © | consommation 3 

Machine à laver importante 


Fig 4 : Légende des plans électriques 1 


Machine à laver 


en = Du — 


tr — / UP 
el — A2 \GN 
L_ 


Réfrigérateur 2 


+ 


Fig 5 : Schéma électrique simplifié de l'étage 1 


Fig 6 : Schéma électrique simplifié de l'étage 2 


Sur les douze interrupteurs différentiels du tableau électrique donné ci-dessous, 
deux ne desservent aucune pièce, comme observé sur le schéma unifilaire 
donné en annexe 1. 


Fig 7 : Photo du tableau 
électrique 


1) Elements de protection 


A) Régimes de neutre, norme NF C 15-100 et rappels de 
dangerosité 


L'installation est en régime de neutre “TT”, c'est-à dire que le neutre et les 
équipements sont reliés à la terre tel qu'illustré ci-dessous. 


Régime TT 


Fig 8 : Régime de neutre TT 


Cependant, contrairement à ce qui est indiqué par le schéma (qui s'applique aux 
gros consommateurs électriques), la maison est en monophasé, c'est-à-dire 
qu'elle ne bénéficie que d'une seule phase. 


La norme NF C 15-100 réglemente les installations électriques en basse tension 
en France. Elle s'applique en revanche uniquement aux installations neuves ou 
entièrement rénovées. Le logement actuel n'a pas été rénové, la norme ne 
s'applique donc pas dessus. En revanche, cette norme doit tout de même être 
respectée au mieux, car elle tient compte non seulement de la protection de 
l'installation électrique et des personnes mais aussi du confort de gestion. Celle- 
ci inclue beaucoup de réglementations, sur la protection des personnes, des 
circuits, contre la foudre, sur les prises de courants, sur les lumières, le tableau 
électrique, la salle de bain et le réseau multimédia. Comme il serait trop long de 
détailler tous les éléments de chacune de ces huit rubriques que contient la 


norme, seuls les éléments principaux de la norme vont être explicités dans les 
parties suivantes. 


B) Protection des personnes 


La norme stipule qu'il faut au moins deux disjoncteurs différentiels de 30 mA 
dans le logement et que maximum huit circuits peuvent être supportés par 
chacun des interrupteurs différentiels. Le disjoncteur différentiel mesure une 
différence de potentiel entre l'amont et l'aval. S'il en détecte une, le courant est 
coupé. || protège les utilisateurs en se déclenchant à partir d'un certain courant 
de fuite, de l'ordre de quelques dizaines de milliampères généralement, car le 
seuil de non lâcher est de l'ordre de 10mA. Il y a effectivement au moins deux 
disjoncteurs différentiels dans le logement (ici douze), mais l'interrupteur 
différentiel orange supporte douze circuits, donc la norme n'est pas respectée. 
Concernant les surintensités, la norme impose l'une des configurations 
suivantes. 


> LE CALCUL DE L'INTENSITÉ 


l 


Intensité Intensité 
interrupteur > disjoncteur 
différentiel d'abonné 


Intensité 
interrupteur 
différentiel 


somme des intensités des 
disjoncteurs pour les circuits 
convecteurs et chauffe-eau 


somme des intensités des 


0 Bb x disjoncteurs pour Les circuits 
; prises, lumières et autres usages 


Fig 9 : Gestion des surintensités imposée par la norme NF C 15-100 
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Dans le cas de la maison, la règle amont est vérifiée, car 
le disjoncteur d'abonné, en photo ci-contre, possède 
une intensité limite de 45 À, supérieure aux 40 A 
annoncés sur l'interrupteur différentiel. 


Fig 10 : Disjoncteur d'abonné 


Concernant la prise de terre, la norme préconise l'un 
des raccordements suivants. 


SOLUTION 1 SOLUTION 2 
BOUCLE EN FOND DE FOUILLE PIQUET DE TERRE 


Cuivre nu 
= de25mm° 


Boucle en cuivre nu 
de 25 mm? 


Fig 11 : Raccordements à la terre préconisés par la 
norme 


C'est la solution 2 qui est privilégiée dans le cas de l'habitation, comme observé 
ci-dessous. 


F Sp U VA fl LR AUOR SUR 2 ae 
Fig 12 : Raccordement à la terre de 
l'habitation 
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C) Protection des biens et des circuits 


Les conducteurs électriques rencontrés sont de trois couleurs différentes : rouge 


pour la phase, bleu pour le neutre et :111-/vert pour la terre (ou masse) tel 


qu'illustré ci-dessous. 


Terre (E) 


Neutre (N) BE — 


Phase (L) 


Jaune / Vert 


Bleu 


Rouge 


Marron et autres 


Fig 13 : Identification des conducteurs électriques 


La maison étant assez ancienne, certains fils électriques visibles pouvaient être à 


nus directement en tirant sur l'ampoule, comme on 
l'observe ci-dessous. 

On observe également le fil de terre simplement enroulé 
autour de la phase et du neutre, ceci n'est pas très 


“propre”. Même si cette installation n'est pas en opposition 


avec la norme, il conviendrait a minima de bien placer 
l'ampoule au plafond de manière à ne pas voir les fils. 
Cette mesure permettrait à l'utilisateur de ne jamais 
toucher un fil sous tension en souhaitant par exemple 
changer une ampoule. 


Sur le tableau suivant, on observe les différentes sections 
de câbles requises selon les cas d'usage. On retiendra les 
sections de câbles du tableau, qui donnent des valeurs 
minorantes, car les sections de câbles réelles sont 
difficilement accessibles dans la maison. Les valeurs 
données dans le schéma unifilaire sont donc des valeurs 
théoriques. 


Fig 14 : Une ampoule dont la 


masse est simplement enroulé 
autour des fils 


On notera cependant que certains câbles pour l'éclairage étant à découvert, leur 


diamètre a pu être mesuré, hors tension bien évidemment. Une valeur de 
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d = 3mm a été trouvée. 


2 3 
On en déduit s = "= 2 77 mm? > 1,5 mm? pour l'éclairage. Cette valeur 


d 

4 4 
peut cependant être faussée, car le câble a été mesuré avec gaine isolante. 
Connaissant la vraie section des câbles, il faudrait vérifier si la chute de tension 
entre l'origine de l'installation et tout point d'utilisation est inférieure à 3% pour 
les éclairages et 5% pour les autres récepteurs, condition exigée par la norme. 


La norme impose également différents types de conducteurs suivant les circuits 
protégés. 


INTENSITÉ 
MAXI 
DISJONCTEURS 


SECTION 
MINI FILS 


CIRCUIT 
PROTÉGÉ 


Circuits lumières 


8 points lumineux maxi 
[ES k par circuit 
ps « Lumières 1,5 mm? 16A 


Au moins 2 circuits 
par logement {> 35 m?) 


Circuits prises de courant 
Prises de courant 
avec terre 


Prises de courant 
avec terre dédiées 
cuisine!) 


Circuits spécialisés 


8 prises maxi 
1 6 A par circuit 
: 12 prises maxi 
2,5 mm | 20 A [par circuit 


2,5 mm° 20 A 6 prises maxi 


æ 


FS Volets roulants 1,5 mm? 1 6 A ché a p pour l'ensemble 
= Chauffage ; 1 circuit dédié par tranche 
[| électrique 2,5 mm 20 A de 4500 W 
Lave-vaissell 
S lave-linge, ” | 1 circuit dédié par appareil 
@ | sèche-linge, 2.5 mm° 20 À  l|aumoins3 circuits 
four électrique, par logement 
congélateur... 
©.®| Plaques / :. 
de cuisson 6 mm 32 A 1 circuit dédié 


Fig 15 : Choix des conducteurs en fonction du circuit protégé 


D'après le schéma unifilaire de l'installation, donné en annexe 1, le nombre de 
lampes maximum supporté par un circuit est de cinq et il y a trois circuits 
lumières, la première ligne est donc respectée. Le nombre maximum de prises 
de courants est de 7 (< 8), pour l'étage 2, donc la deuxième ligne est également 
respectée. On continue ainsi pour chaque ligne du tableau. 
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Après vérification sur le Web, toutes les prises de la maison 
sont de type E (photo ci-dessous), donc reliées à la terre. 
Ceci assure une protection optimale des appareils et des 
personnes. 


Toutes les conditions exigées par la norme n'ont bien sûr 
pas été vérifiées, car il faudrait plusieurs dizaines de pages  TypeE 


(protection contre la foudre, dimensions du tableau E, RARE 


Pologne, Slovaquie & 


électriques, hauteurs des prises, lumières...) pour tout 


vérifier rigoureusement. On peut cependant établir qu'en suc ere 
l'état actuel, la maison ne respecte pas la norme NF C 15- f 220-240 V 
100. 


Fig 16 : Caractéristiques de la 
Ceci n'est pas gênant dans la mesure où, comme rappelé au 


prise E 


début, la norme ne s'applique que pour des logements 
neufs ou récemment rénovés (pas le cas de la maison). 


11) Analyse de la consommation électrique 
A) Démarche 


Toutes les valeurs de tensions et courants donnés seront des valeurs efficaces 
uniquement. 

Le circuit électrique est un montage parallèle de 33 dipôles. 

On considère la tension aux bornes de chaque dipôle de V = 230V + 70V àla 
fréquence de f = 50 Hz + 0,05 Hz. 


Pour chaque dipôle récepteur, diverses informations ont pu être recueillies, par 
trois méthodes différentes : les informations constructrices, la mesure directe 
par Wattmètre et la mesure directe par compteur Linky. 

Il faut noter que pour certains appareils, comme les ordinateurs portables, les 
valeurs sont approchées, car le temps de recharge de la batterie était inconnu. 
Dans ce cas précis, le calcul se base sur une approximation du temps de 
recharge journalier obtenue après discussion avec les utilisateurs (mes 
colocataires) et sur une valeur arbitraire du temps de recharge de la batterie de 
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1h. Prenons un exemple sous ces hypothèses : si le temps d'utilisation est de 
3h, cela signifie que trois cycles de batteries de 1h chacun ont été effectués. 


Le temps de fonctionnement journalier des appareils a été approximé en 
considérant le dernier mois (Mai) seulement. Ce temps de fonctionnement a été 
établi par expérience personnelle et par discussion avec mes colocataires. Pour 
certains appareils, comme la hotte aspirante, que nous n'utilisons que très 
rarement, un temps “symbolique” de 0,05h par jour (3 min) a été appliqué. 


Il faut noter que tous les résultats numériques des puissances actives, réactives 
et apparentes, ainsi que les résultats de l'énergie consommée ont été arrondis à 
l'unité afin d'être cohérent vis à vis des incertitudes de mesures. En effet, ces 
incertitudes relatives sont de l'ordre de quelques pourcents et portent le plus 
souvent sur le chiffre des unités. Elles n'ont pas été calculées par souci de 
simplicité. 


1) Les informations constructeurs 


La puissance active de toute les lampes a pu être 
déterminée par lecture des informations constructeur. 
Par exemple, voici en photo la lampe LED de la chambre 


3 ci-contre. Re 2852 


On peut y lire : Pign = 6W 


Fig 17 : Analyse des informations 


d'une ampoule LED 


En guise d'autre exemple, pour le lave-linge, la notice 
du constructeur donnait les informations suivantes. 
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Fig 18 : Analyse de la notice du lave-linge 


Comme cette notice annonce plusieurs valeurs différentes de puissances actives 
consommées par cycle, on prend le programme le plus utilisé dans mon cas et 
celui de mes colocataires. Il s'agit du programme “Eco 40-60”, que l'on considère 
en capacité nominale. On en déduit la puissance active de chaque cycle. 


E 1,08 
cycle = ah 
7 St; 


On va supposer que 2,5 cycles sont effectués par semaine, car nous sommes 
cinq colocataires et nous effectuons environ un cycle chacun toutes les deux 
semaines. Ceci fait une utilisation journalière Uj de : 


5 LD 
U; = nbcycle/semaine *DUréecycle _ 25x (Gr) 


ND jours/semaine 7 


A titre de vérification, le Wattmètre a été branchée sur la machine à laver 
pendant deux semaines et il a affiché la valeur cumulée d'énergie consommée 
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de 5,58 kWh. En faisant le calcul avec les valeurs estimées précédentes, on 
trouve : 


E> semaines = 14X U; X Ecycie = 14 X 1,22 X 0,316 =54kWh 
L'approximation est donc plutôt cohérente avec la valeur expérimentale trouvée. 


D'autres appareils dont la puissance active ne pouvait être calculée, comme la 
hotte aspirante, ont fait l'objet d'une recherche d'information sur le document 
du projet déposé sur Moodle. 


2) Mesure directe par Wattmètre 


Grâce à un Wattmètre numérique (en photo ci-dessous), 
comme vous l'avez compris, que je pouvais brancher sur 
chaque prise, le facteur de puissance, la puissance active, 
la tension et l'intensité de l'appareil utilisé pouvaient être 
déterminés d'une seule fois. 

Cet appareil présentait l'avantage d'être assez précis. En 
effet À Wattmètre = +2% 

Il s'agit de la méthode la plus efficace car elle donne toutes 
les informations en une mesure et elle sera aussi 
considérée comme la plus fiable, car elle s'établit 
directement (sans établir de différence comme la méthode 
2). 


Ainsi, lorsque le choix se présentera entre deux méthodes, 


seule la mesure expérimentale par Wattmètre sera Fig 19 : Wattmètre utilisé 
pour les mesures 


gardée. 


3) Mesure indirecte par compteur Linky 


Grâce à ce compteur, qui rend compte uniquement de la puissance apparente 
soutirée dans la maison, il était possible de déduire la puissance apparente d'un 
appareil en réalisant une différence de lecture. Par exemple, voici à gauche le 


17 


compteur Linky avant l'allumage du 
four et à droite le compteur après 
allumage du four. 

On déduit immédiatement S;x = 
1350 — 208 = 1,142 kVA 

Cette méthode, bien que très rapide, 
ne donne à elle seule pas toutes les 
informations sur le dipôle étudiée, 


comme la puissance active ou réactive. 
Il faut donc la coupler avec la méthode Fig 20 : Détermination de la puissance apparente 


1 afin d'obtenir toutes les d'un four grâce au compteur Linky 


informations. 


B) Méthode de Boucherot 


1) Fonctionnement standard de l'installation 


La méthode de Boucherot a été appliquée directement depuis le classeur Google 
Sheets où toutes les valeurs utiles au calcul ont été renseignées. Le nombre 
d'appareils est de 33 au total. Le tableau complet pièce par pièce donné en annexe 
2 étant très grand, un tableau synthétique, par type d'appareil, est résumé ci- 
dessous, en réalisant des “groupes” d'appareils et en considérant uniquement les 
grandeurs “utiles" aux calculs. La consommation est étudiée pour une journée 
puis sur un mois. 


Energie 
Groupes 2 Utilisation , 
; FP |Q(VAR)IS (VA)! . consommée 
d'appareils jour (h) , 

(W.h.jour-1) 

-24 51 

0 100 

20 80 

155 781 


lampes 
0,882 
fluocompactes 
ECLAIRAGE lampes ; 
incandescences 
lampes LED 0,958 


Aspirateur 


INFORMATIQUE 
Ordinateurs -293 377 
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Box WIFI 


Imprimante 


micro-ondes 


Frigidaire 


GROS 
ELECTROMENAGER 


Four 
Machine à laver 


Radiateur élec 2000 | 8,696 


Fig 21 : Tableau synthétique des appareils électroménagers de l'installation 


La méthode de calcul va désormais être détaillée. 

Suivant chacune des trois méthodes de détermination données dans la partie 
précédentes, des grandeurs ont pu être obtenues pour chaque appareil et 
renseignées dans le tableau, que ce soit la puissance apparente (S), la puissance 
active (P) ou le facteur de puissance (FP). À partir de deux de ces grandeurs, il était 
possible de déterminer les autres, en sachant au préalable si le dipôle était 
inductif ou capacitif. Les relations utilisées, dans un sens ou dans l'autre, sont les 
suivantes. On notera U; le temps en heure d'utilisation journalier de l'appareil et 
E; l'énergie consommée pour une journée. 


Q 


S=VP+Q = = Vs XI 
Ÿ FP  sin(acos(FP)) ÉMCer 


E; = PXxU; 


A titre d'exemple, pour un aspirateur, on obtient grâce au Wattmètre : 


= 278 


Paspi = 7/65 W et F Paspi = 0,98 = Saspi — re . 0.98 
> Qaspi = S X sin(acos(FP)) = 781 x sin(acos(0,98)) 


= +155 VAR (inductif) 
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S 781 
eff aspi y 230 A 


> | 
Finalement on considère que l'aspirateur est utilisé une vingtaine de minutes par 
semaine, soit environ 3 minutes (0,05 h) par jour : U; = 0,05 h. jour”! 
On en déduit immédiatement : E; = P x U; = 765 x 0,05 = 38W.h 


La méthode de Boucherot permet de déterminer les puissances totales, délivrées 
par le réseau (indicées par “r”), lorsque tous les appareils fonctionnent ensemble, 
par les formules suivantes : 


P,=Y".P. 
” S, = 4 FF + QË 
Q =>; Q; 


On obtient une puissance active du réseau totale de P. = 9,25 kW, une puissance 
réactive du réseau totale de Q, = 0,752 kVAR et une puissance apparente totale 
de S, = 9,28 kKVA > Sontrat = 6kVA. Ainsi, si tous les appareils fonctionnent en 
même temps, des frais supplémentaires vont être causés par cette 
consommation hors contrat. 

Pr _ 9,25 


= = 0.997 > FPmin = 0,93 


Ce facteur de puissance est donc acceptable pour EDF. 


Le facteur de puissance total est de : FP. = 


Toujours si tous ces appareils fonctionnaient simultanément, le courant total 


; S. 9280 
serait de : 1, ef Var — 230 — 41,7 A > Lüisjonction = 40 À 


En effet, le disjoncteur différentiel général est de type “C40”, cela signifie que le 
courant maximum qu'il peut supporter de manière continue sans se déclencher 
est de 40 A. Comme cette valeur est dépassée, il risque d'y avoir une coupure de 


courant. 

L'énergie totale pour un fonctionnement standard consommée en kWh 
journalièrement E;;; est obtenue en sommant les consommations de chaque 
appareil. On obtient, pour une journée standard environ : 


Ets; — X;=1E; = 4,29 kWh 


En considérant un prix d'abonnement journalier P, fixe de 0,38 euros pour le mois 
de Mai (obtenue par un suivi de consommation en ligne) et un prix de 
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consommation du kWh P. de 0,2 euros, on en déduit le coût total standard 
journalier C;.; de l'installation, puis le coût total standard mensuel Ci: pour le mois 
de Mai. 


Crsj = Pa + Ets; X P. = 0,38+ 4,29 x 0,2 = 1,24€ .jour | 
Can ie ANS = Oo) 


Ces coûts sont bien sûr des coûts moyens pour un fonctionnement standard de 
l'installation. 


Grâce à excel, le diagramme en camembert de l'énergie électrique consommé par 
catégorie d'appareil est représenté ci-dessous. 


Répartition de la consommation énergétique 
standard jounalière par type d'appareils 


= Eclairage 
= Informatique 


= Petit électroménager 


= Frigidaire “ 
= Gros électroménager hors 
frigidaire 


Fig 22: Diagramme camembert de la répartition énergétique journalière consommée 


par catégorie d'appareils en fonctionnement standard 


On s'aperçoit de la part importante qu'occupe le gros électroménager tout 
compris, à savoir 45% + 25% = 70% de toute l'énergie consommée chaque jour. 
Le frigidaire/congélateur occupe 45% de toute l'énergie consommée, ce qui en fait 
le plus gros consommateur d'énergie électrique de la maison. 


Finalement, le coût de revient mensuel spécifique de chacun des i appareils pour 
un fonctionnement standard Ci;sn peut être calculé comme suit pour le mois de 
Mai : 

Cr 


Par exemple, pour la bouilloire, on obtient : 
Es; — 0,342 kW. h.jour”! a Cssm bouilloire — 0,2 X 0,342 K51—= 2,12 € 


L'ensemble des coûts calculés pour chacun des appareils est résumé dans le 
tableau suivant. 


Coût mensuel (£.mois'- 


Groupes d'appareils 


lampes fluocompactes 

ECLAIRAGE lampes incandescences 
lampes LED 
Aspirateur 


Ordinateurs 
INFORMATIQUE 
Box WIFI 


Imprimante 


micro-ondes 
Frigidaire 
GROS ELECTROMENAGER Four 
Machine à laver 


Radiateur élec 


Fig 23 : Tableau synthétique du coût mensuel de chaque appareil par catégorie en 
fonctionnement standard de l'installation 
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Bien sûr, le total du tableau ne représente pas le vrai coût total car il ne tient pas 
compte du coût de l'abonnement. Pour plus de confort visuel, voici la répartition 
des coûts par type d'appareil sous forme de diagramme en barre. 


Coût mensuel (€.mois'-1) 


Radiateur élec == 
Machine à laver 
Four ms 
Frigidaire 
micro-ondes = 
Hotte aspirante 
Bouilloire = 
Grille Pain = 
Imprimante 
Box WIFI mms 
a 
[] 
pe 
LL] 
= 


Ordinateurs 

Aspirateur 

lampes LED 

lampes incandescences 
lampes fluocompactes 


0 1 23 4 5 6 7 8 9 10 11 12 13 


Fig 24 : Coût mensuel de l'installation par appareil en 
fonctionnement standard de l'installation 


Après le frigidaire-congélateur, la machine à laver et la bouilloire sont donc les 
plus gros consommateurs d'énergie. 


2) Fonctionnement minimal de l'installation 


Pour un type de fonctionnement minimal, on considère que seul le frigidaire 
associé au congélateur est en marche. On a alors, toujours d'après le tableau ci- 
dessus : 

mb = A our 


> Cnj = Ps + Etmj X P: = 0,38+ 1,92 x 0,2= 076€.jour | 
> CE EE ES 


On constate que le prix journalier et mensuel en consommation minimal est assez 

élevé par rapport à la consommation standard. C'est tout simplement parce que 

le frigidaire est l'appareil représentant la plus grosse consommation d'énergie 

électrique du logement. Comme trouvé précédemment, la part d'énergie 
1,92 


consommée par celui-ci est de 1 — 100! 
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La comparaison entre appareils permanents et appareils amovibles n'a pas été 
réalisée car seul le frigidaire-congélateur fonctionne en permanence, avec le WIFI, 
mais celui-ci représente une part très négligeable de la puissance totale (quelques 
Watts). 


C) Méthode du graphe linéaire orienté et Tellegen 


1. Simulation du circuit 


Afin de modéliser le circuit, une extension? Google Doc a été utilisée. Celle-ci 
permet non seulement de réaliser le schéma électrique simplifié de l'installation, 
mais également de le simuler. À partir des résultats obtenus par la méthode de 
Boucherot précédemment, les valeurs des résistances et réactances (définies par 
des bobines ou condensateurs sur le circuit) ont directement été rentrées. Afin 
de vous permettre de réaliser la simulation, le fichier texte de la simulation, à 
mettre dans un bloc note, est donné en annexe 3. Il suffit ensuite de se rendre sur 
la page https://thumbsdb.herokuapp.com/circuit/ puis dans l'onglet «File » 


d'ouvrir le fichier note. Le circuit se charge alors tout seul. 


Le D Le L L L L L L L L L Le L 
1.64 A(mms)76.7 mA(rms} 49.68 mA(rm<)3 47 A(ims) 718.67 mA(rmès) 55222 mA(rmis) 1.43 Affms) 


220.65 mAfrms 3.95 mAËr 
39.98 Alrms) TA 42 EYES LL) A(rms) 8.84 A(ims) 6.77 A(rmms) 5 A(rms) 8.84 A(rms) 


915 


6.SuF 16.6yF 244 S5mH= 29 2yF= 1 2uF 750nF 
L2 L2 L2 L2 L2 


fuocompactes incandescences LEDS aspirateur ordinateurs WIFI 


mprimante 


Le circuit électrique a été tracé pour les quinze appareils du tableau simplifié 
(Figure 21) et une tension de réseau de 230V. 

Le sens de courant est défini selon le sens des flèches grises, qui correspondent 
à des ampèremètres sur le circuit. Les résistances et valeurs des capacités et 
inductances ont été définies à partir du tableau figure 28 (ci-dessous). Les 
valeurs des résistances correspondent simplement à la lecture graphique du 
tableau, tandis que les valeurs des capacités et inductances définies sur le 
schéma électrique s'obtiennent de la façon suivante : 


2 Ce module est disponible ici : https://thumbsdb.herokuapp.com/circuit/ normalement (il faut 
sinon aller dans les extensions de Google Doc et le télécharger gratuitement) en accès libre. Il 
présente de nombreux avantages pédagogiques et est très simple à utiliser. 
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il 1 
C = —— — 
wX  2rfX 
Re 
w  2rf 


Avec «X» la valeur de la réactance définit dans le tableau (fig 28), « f » la 
fréquence du réseau en Hz. Une simple simulation de la tension du réseau en 
fonction du courant total a pu être effectuée ci-dessous. 


Fig 26 : Courant (en jaune) et tension (en gris) en fonction du temps pour le tableau 


simplifié (Fig 21) de l'installation. 


On trouve un courant total de 39,98 A (en haut à droite Fig 26), proche du courant 
total obtenu par la méthode de Boucherot, mais inférieur. Ceci peut être expliqué 
par une tension mesurée à 238 V environ (voir figure 19) et non 230 V (valeur de 
simulation) dans l'habitation. Ainsi, en mettant une tension plus élevée dans la 
simulation, on retrouverait le même résultat, aux incertitudes de calcul prêt. On 
voit aussi sur le graphique le faible déphasage courant/tension, témoignant d'un 
facteur de puissance très proche de 1, conformément à celui trouvé par la 
méthode de Boucherot (0,997). De plus, toutes les valeurs des courants efficaces 
définies Figure 25 sont bien conformes avec ce qui a été obtenue par la méthode 
de Boucherot. Par exemple, prenons le four, un courant efficace de 4,97 À y a été 
trouvé (tableau fig. 21) et « expérimentalement », nous obtenons 5A. Il faut 
cependant noter qu'il est normal de retrouver les mêmes valeurs, car les valeurs 
des impédances entrées dans la simulation ont été obtenues en connaissance de 
cause du courant y circulant (...). 
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Il se trouve que ce module se base, pour le calcul des courants, sur la méthode de 
calcul automatisée offerte par la topologie de graphes linéaires orientés. 
Maintenant que nous avons introduit le circuit, passons à l'application de cette 
méthode des graphes linéaires orientés. 


2. Application de la méthode 


Le calcul via Python a été effectué pour les 33 appareils. Cependant, afin de ne 
pas surcharger la page, celui-ci a été tracé, comme pour la simulation, pour le 
tableau simplifié (Fig 21) des quinze appareils appartenant aux catégories 
Éclairage, Informatique, Gros électroménager et Petit Électroménager. On aboutit 
au graphe linéaire orienté suivant. On définit le potentiel du nœud 1 
arbitrairement à 0. 


Fig 27 : Graphe linéaire orienté de l'installation simplifiée 


Ce graphe est simple, dans le sens où il présente que deux nœuds, tous les 
appareils étant branchés en parallèle sur le réseau (afin de garantir l'autonomie 
des autres en cas de panne). Il y a une arête par groupe de dipôle du tableau Fig 
21 considéré, soit seize arêtes. On déduit les paramètres du graphe : 
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Rang: p= n -1—=2-1-=1 
Nullité: v=n,-(n -1)=16-1=15 


Il est bien plus avantageux d'utiliser la méthode des potentiels aux nœuds, 
sachant qu'il ÿ a quinze fois plus de branches que de nœuds. Ainsi, utiliser la 
méthode des potentiels aux nœuds nécessite de résoudre qu'une seule équation 
à une inconnue alors que la méthode des courants de mailles nécessite de 
résoudre quatorze équations à 14 inconnues. En revanche, une branche 
comporte une source de tension seule : afin d'utiliser tout de même la méthode 
des potentiels au nœuds, il sera nécessaire plus tard d'utiliser l’’astuce” de la 
petite résistance en série de la source de tension. 

La matrice d'incidence étendue est: 


4 = |? 1 1 1 1 L 1 1 1 1 1 1 1 1 1 1 
$ [1 1 1 1 1 1 1 1 1 L 1 1 1 1 1 J 


La matrice d'incidence s'obtient en définissant le potentiel du noeud 1 à 0, on 


retire donc la première ligne : 


A=[-111111111 11 1 1 1 11] 


La matrice des mailles étendue est donnée ci-dessous : 


1-17 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

0 —1L1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 

0 © =k. À 0 0 0 0 0 0 0 0 0 0 0 0 

0 0 0 I 1! 0 0 0 0 0 0 0 0 0 0 0 

0 0 0 0 —I1 1 0 0 0 0 0 0 0 0 0 0 

D 0 0 0 D SEX 0 0 0 0 0 0 0 0 0 

O0 0 0 0 0 0 —1 1 0 0 0 0 0 0 0 0 

ne — 0 0 0 0 0 0 0 I 1 0 0 0 0 0 0 0 
: 0 0 0 0 0 0 0 O0 1 ! 0 0 0 0 0 0 
0 0 0 0 0 0 0 0 0 I I 0 0 0 0 0 

0 0 0 0 0 0 0 0 0 0 I 1! ES 

0 0 0 0 0 0 0 0 0 O0 O0 —-1 1 0 0 0 

0 0 0 0 0 0 0 0 0 0 0 0 —-1 1 0 0 

0 0 0 0 0 0 0 0 0 0 0 0 0 -I 1! 0 

0 0 0 0 0 0 0 0 0 0 0 0 0 0 I I 

1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -1 


On retire la dernière ligne pour obtenir la matrice des mailles. 


27 


Lo =ts 20% 0 0: 0020 ©: 6 EE 0: 0! CE 0 
O S<E, LL 0 01 0: CO À, 0 Où 0 0. 0 0° 6 
DS 0, ET LE 26 nr, 0 0 CE QD C0 CA 0: 0 | 6 "0 
0 0 © =1° À 0. CO 0, 0 0 0 0 D 0 6 0 
00 0 0 0 tb SE 0 OO 0 O0 0 De À 0 EE EH 
Où CO, CO 20: 0. TE SE 00 C0. 0. D “02: 0 © EE -0 
D: 0 006 C6, C0 0, 2 EE C0 (0 0 LE °0: 0 LE 6 
M=10 0 0 0 0 0 0 -1 1 0 O0 0 0 0 O0 0 
0 0 0 0 0 0 0 0 -1 1 0 0 0 0 0 0 
JO 20: 0 6,0 0 0 0 0 TE TL: 0 20 “0: CE LP 
D, c0 OC: 70: A C0 070 0 O0 SEE 2% CE Xe CS 
D 0-02 (0 OÙ -<6 0: 0, 0 Or (DE SE CP 20: 0 
OÙ De 20 CO A0 07 C0 Er Or 0 00 0er =" AE CS 
OS Où. OC CO Æ 07 0 00: C0: 0 ® EN" 0: EL 1: -0 
0: “©@, 0, 0, 2 © 07 “0 0 D V0 D 0 Pet. 


On définit ensuite pour chacune des quinze arêtes un composant générique 
complexe comme ci-contre 


Ve 


Hormis la tension E (-230V) du réseau 
présente seule dans une branche, il n'y a pas 
de sources de courants dans les branches. Le 
vecteur Is est donc constitué d'une colonne 


avec seize lignes remplies de 0. Par soucis de 


place, ces vecteurs sont représentés en ligne complexe 
ci-dessous. 
Vs=[E 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0] 
Is=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 


Comme il y a une tension seule dans une branche, cela va conduire l'admittance 
Y à exploser dans la résolution du système linéaire. Pour cette raison, on introduit 
arbitrairement une petite impédance R = 0,019 en série avec la source de 
tension. Toujours grâce au tableau simplifié, les valeurs des impédances 7, et 
admittances Y, ont pu être calculées. 
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Groupes Résistance | Réactance , Admittance complexe 
: | Impédance complexe (Q) 
d'appareils (Q) 


0,000850661625708885 
lampes 915,22491349481- 
-488 ; +0,00058273193685343 
fluocompactes 488,119953863899i > 
lampes 
ÉCRAN | P 529 0,00189035916824197 
incandescences 
634,178256611166- 0,0014461247637051+0, 
lampes LED -191 : : 
190,667972575906i 00047407901712641i 


70,2479338842975+14,264 | 0,01367157647/705882- 


Aspirateur - , 
4513668515 0,00288825860878766i 


88,2796546916588- 0,00448015122873346+ 
109,138982382515i 0,013997026960638i 


Ordinateurs 


INFORMATIQUE More 1433,90117647059- 0,000160680529300567 
OX 
2620,65899373516i +0,0012745921111484i 


0,000086956521739130 
Imprimante 1840-4215,96963935937i | 4+0,0012452636632516 
3i 


Grille Pain 66,3212758652976 0,0150781176470588 


PETIT Bouilloire 25,8048780487805 0,0387523629489603 


ELECTROMENAGER 
272,515151515152+168,88 | 0,00265122873345936- 


Hotte aspirante 9722280579 0,00227404221885373i 


30,7969897358919+13,669 | 0,0271266540642722- 
2752785943 0,0143516295807253i 


micro-ondes 


262,450125+323,51982480 | 0,00151228733459357- 
79 0,00469658418372534i 


Frigidaire 


GROS 
ELECTROMENAGER Four 46,3222416812609 0,0215879017013233 


154,280506329114+44,998 | 0,00597353497164462- 


Machine à laver ; ; 
4810126582i 0,00189017826898088i 


Radiateur élec 26,45 0,0378071833648393 


0,000070485838358548 
6395,77310332802- 4+0,0001725719456547 
7059,21378717027i A6i 


Fig 28 : Tableau des impédances et admittances complexes par appareils 
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Les résistances et réactances ont été calculées comme suit : 


Fe 


Z=R+jX 
X = 


lo l* 


La puissance complexe a été calculée par la fonction COMPLEX(Re:; Im) de Google 
Sheet prenant en paramètre la partie Réelle, R et la partie imaginaire, X. 
Finalement, comme Google Sheet ne permet pas d'inverser un nombre complexe 
directement, la formule suivante à été utilisée : 


Es 2 
_R+jX  R2+X2 


À partir de ces données, la matrice d'admittance de dimensions 75 x 15 a été 
calculée. 


Passons désormais au calcul des courants et puissances dans les différents 
dipôles grâce à la méthode des potentiels aux nœuds. Afin d'appliquer cette 
méthode malgré la présence d'une source de tension seule dans une branche, 
une petite résistance de 0,01 N a été mise en série de la source de tension. 


Afin de réaliser les calculs, le programme Python déposé sur Moodle a été repris, 
adapté et commenté. Celui-ci est disponible en annexe 4 ainsi que les résultats 


commentés en annexe 5. 


Les résultats obtenus par la méthode des graphes linéaires orientées sont 
cohérents avec ceux obtenus par Boucherot et par la simulation. 
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IV) Comparaison à la consommation réelle et 
perspectives d'amélioration 


À) Comparaison à la consommation réelle 


Enedis et GRDF, en partenariat avec le Ministère de la Transition écologique, 
et l'ADEME se sont associés pour lancer le «Energy Data Challenge »° 
permettant de récompenser les services les plus innovants en termes de 
données de consommation et de production d'énergie fournies aux 
consommateurs. 

Dans ce cadre, deux services de visualisation de données ont été utilisés. 


1) Le service du distributeur Enedis directement 
2) Le service HelloWatt 


Grâce au numéro de Point de livraison d'électricité (PRM) donné sur le compteur 
Linky ou une facture, il a été possible d'avoir les données de consommation de la 
maison en temps réel. 

La consommation suivie pour le mois de Mai est représentée ci-dessous. 


38€ 125 KWh 3,9 kgCO: (@) < Mai 2023 > Année Mois Jour 


01 03 05 07 09 11 13 15 17 19 21 23 25 27 29 31 
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Fig 29 : Consommation électrique de la maison pour le mois de Mai 


On calcule l'énergie moyenne journalière consommée par jour sur le mois de Mai : 


125 
E yai th — 31 — 4,03 kWh 


Cette valeur est à comparer avec la valeur trouvée par la méthode de Boucherot 


3 Ce challenge est détaillé sur cette page : https://energiedatachallenge.fr/ 


cal 


Epai exp = 429 kWh 


L'écart relatif est donc de E, = lEmaich=Emaiexpl _ 6,5 %. 

Emaith 
Cet écart est cohérent avec les incertitudes de mesure, notamment sur les temps 
d'utilisation, qui ont dû être approximés et résultent uniquement de discussions 
avec mes colocataires. Vu le nombre d'appareils de l'installation, on peut dire avoir 
tout de même trouvé une valeur cohérente d'énergie consommée. 


Le prix de l'électricité payé à cinq pour le mois de Mai (38 euros) est plus que 
raisonnable et on peut dire que nous ne sommes pas de très grands 
consommateurs. En effet, la consommation moyenne par colocataire est de 25 
KWh (125/5 car nous sommes 5), or la consommation moyenne d'un français est 
de 185 kWh par mois! 

Finalement, HelloWatt affiche aussi les émissions approximatives de CO2 sur le 
mois de Mai. Elles sont de 3,9 kg environ ce qui correspond à un trajet de 1645 
km en TGV d'après le site. En ramenant ces émissions par personne, on obtient 
3,9/5 soit 0,78kg CO2 par mois par personne, ce qui est très faible ! A titre de 
comparaison, un humain rejette en respirant entre 0,7 et 0,9kg de CO2 chaque 
jour. 

Le contrat souscrit S = 6kVA paraît parfaitement adapté à nos besoins, car cette 
valeur n'a jamais été dépassée d'après Enedis. En effet, la valeur la plus haute 
mesurée serait de 5,2 kVA, valeur atteinte le 28 Février, sans doute à cause du 
radiateur électrique en marche (Hiver). 


1,31€ 4,47 KW 160 gCO; (@) < Mercredi 10 mai 2023 > Année Mois Jour 


Oh 3h 6h 9h 12h 15h 18h 21h Oh 
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Fig 30 : Consommation électrique pour la journée du Mercredi 10 Mai 2023 


En analysant le profil type de consommation sur une journée, on 
obtient naturellement un pic de consommation vers midi, du fait de l'utilisation 
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des appareils électriques en cuisine, comme la bouilloire, le micro-ondes ou le 
four et également en soirée, car non seulement tous les étudiants rentrent de 
cours et branchent leur ordinateur, mais le repas surtout sollicite les mêmes 
appareils que le midi. 


En analysant une journée où nous étions tous en vacances, il est possible de 
trouver la consommation du seul appareil permanent : le frigidaire-congélateur, 
a fortiori aussi le plus gros consommateur d'énergie de la maison. 


On trouve, toujours grâce à HelloWatt, une consommation journalière de 2,4 KWh 
lorsque le frigidaire fonctionne seul. La valeur trouvée dans le tableau 21 était de 
1,92 KW. La valeur a donc été sous-estimée. 


B) Perspectives d'améliorations 


Il aurait été intéressant, avec plus de temps, d'établir les différentes incertitudes 
sur les mesures des puissances et énergies. Ainsi, on aurait pu voir à quel point 
les résultats obtenus en termes de consommation énergétique et de temps 
d'utilisation des appareils sont cohérents avec la réalité. 


Par ailleurs, avoir les données de consommation de l'année précédente (anciens 
colocataires) sur le site d'Enedis aurait permis de les comparer avec les valeurs 
actuelles. 


Le frigidaire-congélateur est l'appareil consommant le plus dans la maison. Il 


serait intéressant donc, à terme, de passer à un appareil très basse 
consommation (label A+). 
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Conclusion 


Ce projet m'a permis d'appliquer les connaissances théoriques à un cas concret, 
me familiarisant ainsi avec les différents aspects de la protection et de la 
distribution dans les réseaux domestiques. Désormais, je sais appliquer des 
méthodes quasi-systématiques pour l'étude des installations domestiques. Les 
méthodes ne doivent pas être appliquées indépendamment. La méthode de 
Boucherot, très pratique et facile à mettre en place, doit être combinée avec la 
méthode automatique des graphes linéaires, afin de bénéficier de leurs avantages 
respectifs. 


De plus, afin de réduire sa consommation énergétique, il est nécessaire de savoir 
où agir. Ce projet donne des clés essentielles pour identifier les appareils les plus 
coûteux et devenir un consommateur responsable. 


Finalement, d'un point de vue sécurité, ce projet m'a permis de comprendre à quel 
point la norme NF C 15-100 pouvait être contraignante et difficile à vérifier pour 
une simple installation domestique, mais à quel point elle était essentielle pour 
garantir la sécurité des biens et des personnes. 


La prochaine étape est d'étudier une installation électrique industrielle, 


fonctionnant généralement en triphasé, qui permettrait d'intégrer pleinement les 
bases de l'électro-technique. 


34 


V) Annexes 


-electrotech 


A1 : Schéma unifilaire de l'installation 


osccsdiccc sectes BUT: 
LL 


20LB2S2pUEIU| 8 adUUE|.: S2UPISSPUBIUI 


Mana tu tend etaenteece cs RE. 


+ 


CREER | PTE I “DNEPELT. - 
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A2 : Bilan global des appareils 


lampe fluocompacte 


lampe fluocompacte 


lampe entrée 


e incandescence 
Entrée 


+ 


Couloirs 
+ 


Salon 


lampe LED entrée 
aspirateur 
lampe LED salon 
lampe LED escalier 


Box WIFI 


micro-ondes 
Frigidaire 
Grille Pain 
Cuisine Bouilloire 


Four 


Hotte aspirante 


lampe LED 


lampe LED 
Salle de bain lampe LED 


Machine à laver 


lampe LED 1 


lampe LED bureau 1 
lampe LED 4 


lampe LED 2 


Radiateur élec 5 
Imprimante 5 
lampe Fluo-compacte 5 
Chambres PCI 
PC2 
PES 
PC 4 


CS 


lampe Incandescence 3 


lampe LED 3 


lampe LED 3 


| TOUS LES APPAREILS 


A3 : Code Simulation‘ 


$ 1 0.0000049999999999999996 10.20027730826997 50 5 43 
v 48 352 48 288 0 1 50 325000.5 

w 48 208 96 208 0 

w 48 352 96 3520 

8 48 352 48 400 0 

r 96 256 96 3040915 

c 96 304 96 352 0 0.0000065 11.360950672489698 
r 160 272 160 352 0 529 

r 208 256 208 304 0 634 

c 208 304 208 352 0 0.0000166 -10.228408305360361 
r 272 272 272 304070 

1272 304 272 352 0 0.0445 3.7188456084032557 
r 336 240 336 304 0 88 

c 336 304 336 352 0 0.000029199999999999998 115.38582856307848 
r 400 240 400 304 0 1434 

c 400 304 400 352 0 0.0000012099999999999998 173.57517748498785 
r 480 240 480 304 0 1840 

c 480 304 480 352 0 7.5e-7 201.91473316830957 
r 560 240 560 352 0 66.3 

r 608 272 608 352 0 26 

r 656 240 656 304 0 273 

w 96 208 160 208 0 

w 160 208 208 208 0 

w 208 208 272 2080 

w 96 352 160 3520 

w 160 352 208 3520 

w 208 352 256 3520 

w256 3522723520 

w 272 208 336 208 0 

w 272 352 336 3520 

w 336 208 400 208 0 

w 400 208 480 208 0 

w 480 208 560 208 0 

w 560 208 608 208 0 

w 608 208 656 208 0 

w 336 352 400 3520 

w 400 352 480 3520 

w 480 352 560 3520 

w 560 352 608 352 0 

w 608 352 656 3520 

r 736 272 736 304 0 31 

w 656 208 736 208 0 

w 656 352 736 3520 

1656 304 656 352 0 0.53 0.578309487034596 
1736 304 736 352 0 0.044 6.410429826634716 
r 816 240 816 288 0 262 

1816 288 816 352 0 1.03 0.20278222560847028 
w 736 352 8163520 

w 736 208 816 208 0 

r 896 272 896 352 0 46 

r 944 240 944 288 0 154 

r 992 256 992 336 0 26 

1944 288 944 352 0 0.14 1.557241157701037 
w 816 208 896 208 0 

w 896 208 944 208 0 

w 944 208 992 208 0 

w 944 352 896 3520 

w 816 352 896 3520 

w 944 352 992 3520 

w 992 336 992 3520 

w 992 224 992 208 0 

x 57 373 121 376 4 10 fluocompactes 

x 128 375 192 378 4 9 incandescences 

x 201 375 230 378 4 12 LEDS 

x 248 373 292 376 4 10 aspirateur 

x311 375 371 378 4 12 ordinateurs 

x 392 379 417 382 4 12 WIFI 

x 451 380 510 383 4 12 imprimante 

x 533 380 578 383 4 10 Grille-pain 

x 591 381 632 384 4 10 Bouilloire 

x 645 379 667 382 4 10 hotte 

x 709 379 764 382 4 10 micro-ondes 

x 804 379 823 382 4 10 frigo 

x 921 377 963 380 4 10 lave-linge 

x 887 376 904 379 4 10 four 

x 977 377 1016 380 4 10 radiateur 

370 48 288 48 20811 

370 992 224 992 25611 

370 944 208 944 240 1 1 

370 896 208 896 27211 

370 816 208 81624011 

370 160 208 16027211 

370 208 208 208 25611 

370 96 208 96 256 11 

370 272208 27227211 

370 336 208 336 24011 

370 400 208 400 240 1 1 

370 480 208 480 240 11 

370 560 208 560 240 11 

370 608 208 608 272 11 

370 656 208 656 24011 

370 736 208 73627211 

x 18 179 1085 182 4 24 0,222\s;\s0,435\s;\50,347\s;\53,3\s;\51,638\s;\s\50,077\s;\50,050\s;\53,58\s;\58,91\s;\50,72\s;\56,83\s;\50,55\5;\54,97\s;\51,43\s;\58,696 
x 8 138 1091 141 4 24 Valeurs\sthéoriques\sdes\scourants\strouvés\sdans\sle\stableau,\sdans\sl'ordre\s(ä\scomparer\savec\sla\ssimulation): 
0 76 8 0 20487 640 12.8 0 2 763 


4 (à copier-coller sur https://thumbsdb.herokuapp.com/circuit/) 
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A4 : Programme Python Utilisé 


import numpy as np #Importation du module de calcul scientifique numpy utile pour les calculs 


d'algèbre linéaire 


texte a afficher = " Mini-projet : Calcul d'une installation électrique domestique " 
#Affichage du titre du programme 
print ( (lén(texte. & affichér) + SjAÆMAeRE MAQUE DRMAN E Léxte à. BÉFICher #.2# 108 À 


+ Mn" + (lén(texte à afficher) + 3)*"4") 


#---—- Paramètres du réseau----- 
E = 230 #11 s'agit de la tension de la source 
#----—- Equations d interconnexion----- 
A = np.ones((1,16))*(1) #Calcul de la matrice d'incidence réduite (partie IIIc) : une matrice 


de taille 16 lignes et 1 colonne est remplie automatiquement 


A[0O,0]=-1 #Le premier coefficient est négatif dû à la présence de la source de tension 


rho = np.linalg.matrix rank(A) #Calcul du rang de la matrice À : nombre de lois des noeuds 


indépendantes (-1 ici) 


print( "\nMatrice d'incidence reduite À = \n", A, "\n de rang rho = ", rho) #Affichage de la 
matrice calculée 
M = np'array([[L, L, 0, 0; 0, 0; 0,0; 0:. 0,0, O0, 0, 0,0; O1, #Entrée de la matrice des 
mailles réduite (partie IIlIc) 

Dieu, 1; 07 07-10, 00,0, 0010 00 D, 0 00101, 

D 0, tr 1 0720, 0% 07 01000, oi 0 +0; 0: “D. 

0.0, 0, ed; 2, 0 04-00, 0, 20100, 0,07; “0, ©: 01, 

0. 0, 0, 07 =, %, 0,.07-00;,-.0,- 0, Q, 0,0, 0. 01, 

0, 0,. 0, 0, 0, 21, 2, 0, 0,0, 0, ©, 0, 0, 0, OT, 

0, 0, 20, 00, 0, 0, =1h , 00, 10; 0, -0, 0, 0; 0, "D, 

D; 0 0,0; 0: 20, 0, 2€, 4, 07 0% 0 G,«0; 0 401, 

0 0, 0,v0, 00, 0, D el 10,0; 0;#0, 0} 01, 

0 s0m 0,10, 60 0, 0,10, St, LE, 0: 0, 10, 00,011 

0, 0, 0: 0, 0, 0: O7 6:00 0, =1, I, 6; "0; 0: 01, 

0, “0; 20, 0, 0: ©: ©, 40; 10, (0, 0 =1, L, 0; 0,01, 

D: 10 0,0% 0 0, 00 0:00; 0-16 û, 1; 1; 0 401. 

D, 05 0"0,. 0 0-0, 0 076107 0e 65e, dd 1407, 

O0, «0, 0, 0, 0, 0, 00, 0, 0%. 0:10, 0, 0, 0, =1,; 111) 
nu = np.linalg.matrix rank (M) #Rang de la matrice des mailles : nombre de lois des 
mailles indépendantes (-15 ici) 
print ("\nMatrice des mailles reduite M : \n", M, "\n de rang nu = ", nu) # Affichage de la 
matrice des mailles réduite 
print ("\nVerification du theoreme de Tellegen : A*M.T = \n", A.dot(M.T)) # Le produit 


scalaire des vecteurs tensions/courants doit être nul, c'est l'équivalent matriciel ici 


Vs = np.zeros((16,1),dtype = complex) #On définit le vecteur des tensions sources comme une 


colonne à 16 lignes remplies de 0. 


Vs[0,0]=-E # Le premier coefficient vaut cependant la tension de la source (non nul) 
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Is = np.zeros((16,1),dtype = complex) # Il n'y à pas de générateur de courant donc on 


définit un vecteur à 16 lignes remplies de 0 


#La matrice d'impédance est définie par un vecteur ci-dessous 

Z_ temp = np.array([0.001,915.22491349481-488.119953863899j,529,634.178256611166- 
190.6679725759063,70.2479338842975+14.2644513668515;j,88.2796546916588- 
109.138982382515;,1433.90117647059-2620.65899373516;j,1840- 
4215.96963935937;j,66.212758652976,25.8048780487805,272.515151515152+168.889722280579;j,30.79698 
97358919+13.6692752785943j,262.450125+323.5198248079j,46.3222416812609,154.280506329114+44.998 
4810126582j,26.45], dtype = complex) 

Z = np.diagflat(Z temp) #Ce vecteur est transformé en matrice diagonale par cette fonction. 

Y temp = [1/r for r in Z temp] #Chaque coefficient du vecteur Z temp est inversé pour obtenir 


la matrice d'admittance. 


Y = np.diagflat(Y temp) #Similairement, le vecteur Y temp est transformé en matrice diagonale 
# —---- Methode des potentiels aux nœuds ----- 

Upsilon = np.linalqg.solve(A.dot(Y).dot(A.T), A.dot(Y).dot (Vs) - A.dot(Is)) #Résolution du 

système linéaire (1X1 ici) pour trouver le vecteur des potentiels aux noeuds 

V = A.T.dot (Upsilon) #On remonte sur les tensions 

T = Y,dot(V) + 15.= Y.dét (VS) #On en déduit les courants par l'équation constitutive 

print ( "\nResultat de la methode des potentiels aux noeuds :\n", "\nUpsilon = \n", #Affichage 


des vecteurs V et I trouvés par cette méthode 
Upsilon, "\n\nv = \n", V, "\n\nI = \n", I) 


# —---- Methode des courants de mailles ----- 
lota = np.linalg.solve(M.dot (Z).dot (M.T), M.dot(Z).dot(Is) - M.dot (Vs)) #Résolution du système 
linéaire (15X15 ici) pour trouver le vecteur des courants de mailles 
TI = M.T.dot(lota) #On remonte sur les courants 
V = Z.dot(l) + Vs - Z.dot(Ils) #On en déduit les tensions par l'équation constitutive 
print ("\nResultat de la methode des courants de mailles :\n", "\nlota = \n", Iota, #Affichage 


des vecteurs I et V trouvés par cette méthode 
MONDTI = Nat, LT ' MMNAV = Ni, VW ) 


# —---- Calcul des pssances et verification du theoreme de Tellegen ----- 
Scomplexe = I.conj() * V # Calcul de la puissance complexe 
print ("\nPuissance complexe = \n'", Scomplexe, 
"\n\nVerification du theoreme de Tellegen : V*I.T =\n", I.conj().T.dot(V) ) # Le produit 


scalaire des vecteurs tensions/courants doit être nul, 


print ( "\nValeurs efficaces des tensions (V) = \n'", np.absolute(V), # La fonction absolute 
renvoie le module de chaque complexe constitutif du vecteur V (-valeur efficace) 
"\nAngle des tensions (degres)= \n", np.angle(V, deg=True) ) # Calcul du déphasage de la 


tension par la fonction angle renvoyant l'argument de chaque complexe de V 


print ( "\nValeurs efficaces des courants (A) = \n'", np.absolute(l), # La fonction absolute 
renvoie le module de chaque complexe constitutif du vecteur I (-valeur efficace) 

"\nAngle des courants (degres)= \n", np.angle(Il, deg=True) ) # Calcul du déphasage de la 
tension par la fonction angle renvoyant l'argument de chaque complexe de I 

cond AYAT = np.linalg.cond(A.dot(Y).dot(A.T)) 

print (cond AYAT) #Affichage du conditionnement : il vaut 1 
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A5 : Résultat du programme Python 


AMI = 


[[O. O. O. 0. O. O. O. O. O. O. O. O. O. 0. O.]] 


Upsilon = 
[[229.96019374+0.00302247j]] 


V= 
[[-229.96019374-0.00302247i] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247j] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247j] 
[229.96019374+0.00302247j] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247j] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 

[229.96019374+0.00302247j]] 


1= 
[[3.98062626e+01-3.0224/7428e+00j] 
[1.95616941e-01+1.04332338e-01;] 
[4.34707361e-01+5.71356197e-06j] 
[3.32549817e-01+9.99870638e-02j] 
[3.14392676e+00-6.38358532e-01;] 
[1.03023970e+00+1.27370611e+00j] 
[3.69492380e-02+6.75321124e-02;j] 
[1.99959364e-02+4.58180887e-02j] 
[3.47304958e+00+4.56479135e-05;] 
[8.91150089e+00+1.17128020e-04;j] 
[6.09682039e-01-3.77835901e-01;j] 
[6.23808/702e+00-2.76868314e+00j] 
[3.47771523e-01-4.28683201e-01;j] 
[4.96435806e+00+6.52488777e-05;] 
[1.37368053e+00-4.00637229e-01;] 
[8.69414721e+00+1.14271239e-04j]] 


Le théorème de Tellegen est bien respecté, c'est à dire que 
la somme algébrique des puissances instantanées dans le 
circuit fermé de l'habitation est nul. 


Le potentiel au nœud 1 est bien de 230V. 


La tension de la première branche est négative car il s'agit 
de la branche où se trouve la source, or onesten 
convention récepteur, d'où le « -». Toutes les tensions ont 
bien une valeur absolue proche de 230, la tension efficace 
du réseau. La partie imaginaire obtenue résulte de légères 
erreurs de calcul. 


Les courants complexes possèdent bien une réactance 
négative et positive respectivement pour des appareils 
capacitifs et inductif. 


Ainsi, les lampes fluocompactes (L2), lampes LED (L4), les 
ordinateurs (L6), la Box WIFI (L7) et l'imprimante (L8) ont 
clairement un comportement capacitif, d'où un courant en 
avance sur la tension, donc une partie imaginaire positive. 


En revanche, l'aspirateur (L5), la hotte aspirante (L11), le 
micro-ondes (L12 : pour cet appareil, c'est une supposition, 
car le comportement varie selon les appareils), le frigidaire 
(L13) ou la machine à laver (L14) ont un comportement 
plutôt inductif, d'où un courant en retard sur la tension, 
donc une partie imaginaire négative. Pour les autres 
appareils, purement résistif, ils présentent généralement 
une partie imaginaire ridiculement faible. 
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lota = 
[[39.80626261-3.02247428e+00j] 
[39.61064567-3.12680662e+00j] 

[39.1759383 -3.12681233e+00j] 
[38.84338849-3.22679940e+00) 
[35.699461 72-2.58844086e+00) 
[34.66922202-3.86214697e+00) 
[34.63227278-3.92967908e+00) 
[34.61227684-3.9754971/7e+00)j 
[31.13922726-3.97554282e+00) 
[22.22772637-3.97565995e+00) 
[21.61804433-3.59782405e+00) 
[15.37995732-8.29140910e-01) 
[15.03218579-4.00457709e-01); 
[10.06782774-4.00522958e-01) 
[ 8.69414721+1.14271239e-04) 


] 
] 
] 
] 
] 
] 
] 
] 


dd NI LU 


] 


1= 

[[3.98062626e+01-3.0224/7428e+00j] 
[1.95616941e-01+1.04332338e-01;] 
[4.34707361e-01+5.71356197e-06j] 
[3.3254981/7e-01+9.99870638e-02j] 
[3.14392676e+00-6.38358532e-01;] 
[1.03023970e+00+1.27370611e+00j] 
[3.69492380e-02+6.75321124e-02j] 
[1.99959364e-02+4.58180887e-02j] 
[3.47304958e+00+4.56479135e-05;] 
[8.91150089e+00+1.17128020e-04;j] 
[6.09682039e-01-3.77835901e-01;j] 
[6.23808/702e+00-2.76868314e+00j] 
[3.47771523e-01-4.28683201e-01;j] 
[4.96435806e+00+6.5248877/6e-05;] 
[1.37368053e+00-4.00637229e-01;] 
[8.69414721e+00+1.14271239e-04j]] 


V= 
[[-229.96019374-0.00302247i] 
[229.96019374+0.00302247) 
[229.96019374+0.00302247) 
[229.96019374+0.00302247) 
[229.96019374+0.00302247) 
[229.96019374+0.00302247) 
[229.96019374+0.00302247) 
[229.96019374+0.00302247) 
[229.96019374+0.00302247) 
[229.96019374+0.00302247) 


] 
] 
] 
] 
] 
] 
] 
] 
] 
[229.96019374+0.00302247)] 


Les courants de mailles ont tous une partie réelle positive, 
cohérente avec la convention récepteur choisie. 


Les résultats sont rigoureusement identiques à ceux 
trouvées par la méthode des potentiels aux nœuds, malgré 
une taille 15 fois plus grande du système à inverser. La 
légère différence qu'il peut y avoir à la dernière décimale 
est liée à la précision de calcul. 


Là encore, les résultats sont identiques et les mêmes 


commentaires s'appliquent pour cette méthode. 
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[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247)] 
[229.96019374+0.00302247j]] 


Puissance complexe = 
[[-9.15384673e+03-6.95169085e+02j] 
[4.49844250e+01-2.39916933e+01;j] 
[9.99653889e+01-2.16840434e-19j] 
[7.64735225e+01-2.29920394e+01;j] 
[7.22976078e+02+1.46806554e+02j] 
[2.36917972e+02-2.92898589e+02;] 
[ 8.49705805e+00-1.55295860e+01;j] 
[4.59840789e+00-1.05362761e+01;j] 
[7.98663155e+02+0.00000000e+00j] 
[ 2.04929047e+03+0.00000000e+00j] 
[ 1.40201458e+02+8.68890598e+01); 
[1.43450333e+03+6.36705765e+02) 
[7.99723111e+01+9.85811231e+01)j 
[1.14160474e+03+0.00000000e+00)j 
[3.15890629e+02+9.21347667e+01) 


] 
] 
] 
] 
] 
[ 1.99930778e+03+0.00000000e+00j]] 


Verification du theoreme de Tellegen : V*I.T = 


[[-5.58962827e-10-1.51287594e-10j]] 


Valeurs efficaces des tensions (V) = 
[[229.96019376] 

[229.96019376] 

[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 
[229.96019376 


UN UN UN UN UN UN UN UN UN UN Lu Lu Eu nu 


] 


Angle des tensions (degres)= 
[[-1.79999247e+02] 
[7.53065203e-04] 

[ 7.53065203e-04] 
[ 7.53065204e-04] 


La puissance complexe trouvée est en accord avec les 
données du tableau figure 21. La puissance apparente 
négative (sans sens physique) de la première ligne 
correspond à celle de la petite résistance en série (0,01 
Ohm), qui voit à ses bornes une tension « négative ». 
(convention récepteur). 


Bien évidemment, la même remarque que précédemment 
est à faire : si l'appareil est à tendance capacitive, sa 
puissance complexe aura une partie imaginaire négative, et 
inversement. 


Le théorème de Tellegen est bien vérifié : la somme des 
puissance instantanées de tous les composants (y compris 
la source) est bien nulle, aux incertitudes de calcul de la 
machine près (de l'ordre de 10-16). 


Les valeurs efficaces correspondent aux modules des 
tensions complexes précédemment trouvées. Elles sont 
donc toutes positives, et proche de la tension efficace du 
réseau (230V). 


Les tensions, définies comme référence de phase, ont un 
déphasage quasi nul. 
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[ 7.53065204e-04 
[ 7.53065204e-04 
[ 7.53065202e-04 
[ 7.53065201e-04 
[ 7.53065202e-04 
[ 7.53065202e-04 
[ 7.53065202e-04 
[ 7.53065202e-04 
[ 7.53065202e-04 
[ 7.53065202e-04 
[ 7.53065202e-04 
[ 7.53065202e-04 


UN UN UN UN UN UN UN UN UN UN nu Lu 


] 


Valeurs efficaces des courants (A) = 


[[39.92084535] 
[0.22170075] 

[0.43470736 
[0.34725609 
[3.20807997 
[1.63820667 
[0.07697943 
[ 0.04999135 
[3.47304958 
[ 8.91150089 
[0.71726714 
[ 6.82490556 
[0.55200935 
[ 4.96435806 
[1.43091173 


] 
] 
] 
] 
] 
] 
] 
] 
] 
] 
] 
] 
] 
[ 8.69414721]] 


Angle des courants (degres)= 


[[-4.34211485e+00] 
[2.80732400e+01] 
[ 7.53065203e-04] 
[1.67343473e+01] 
[-1.14775879e+01] 
[5.10322677e+01] 
[6.13153511e+01] 
[6.64225746e+01] 
[ 7.53065202e-04] 
[ 7.53065202e-04] 
[-3.17875776e+01] 
[-2.39333373e+01] 
[-5.09491244e+01] 
[ 7.53065202e-04] 
[-1.62594516e+01] 
[ 7.53065202e-04]] 


On retrouve presque exactement les mêmes valeurs que 
celles du tableau figure 21 ! Par exemple, pour les lampes 
fluo-compactes (L2), on trouvait lexp = 0,222 A d'après le 
tableau. Ici, on trouve | =0,2217 À, soit des valeurs très 
proches. 


On trouve en revanche un courant total It = 41,7 À, tandis 
que celui circulant dans la petite résistance fictive de 0,01 
Ohm est de 39,92 A. L'écart relatif entre les deux valeurs est 


donc de : 


: ES = 2 en 
= | —— | x — 
: 41,7 oo 


Cet écart s'explique par deux raisons : 

1) La présence de cette légère résistance, qui 
consomme un peu de courant. 

2) La raison principale reste la valeur approximative 
de la tension du réseau, considérée à 230 V alors 
qu'elle était proche de 238 V dans là maison. En 
refaisant les calculs avec 238V, on trouve It=41,3 A, 
beaucoup plus proche de la valeur trouvée par le 
th. De Boucherot. 


Les angles des courants sont respectivement positif s'il 
s'agit d'un appareil capacitif et inductif sinon, 
conformément à ce qui a été dit. On retrouve bien ces 
résultats ici. Par exemple, les lampes fluocompactes (L2) 
ont bien un angle des courants &i positif, le courant est 
donc bien en avance sur la tension dans ce cas. 
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